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ABSTRACT 

Earthquake data from California from the period of 
January 1, 1969, to December 31, 1971, has been analyzed 
by two methods. The first method compares the occurrence 
of earthquakes with the Synodic, Draconic and Anomalistic 
Lunar Periods. Statistical tests to determine if the earth- 
quakes are uniformly distributed over each of the three 
lunar periods are presented. The second method compares 
three components of the tidal force vector and their 
derivatives at the times of earthquake occurrences. 
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I. INTRODUCTION 



A lunar correlation with the time o£ the occurrence 
o£ earthquakes as reported by Allen [Re£. 1] has been one 
o£ the methods used in trying to determine the causes o£ 
earthquakes £or the better part o£ a century. Knopo££ 

[Re£. 3] reported that it may be that earthquakes have a 
greater tendency to be triggered when the sun and the moon 
are nearly aligned with the earth. An interesting observa- 
tion which was reported by Wigand [Re£. 16] was that the 
sun and the moon were almost aligned with the earth at 
the time o£ the San Fernando Valley Earthquake o£ 

February 9, 1971. It is Further contended by Wigand 
[Re£. 16] that there is strong evidence that some earth- 
quakes are related to tidal Force, and some are actually 
triggered by the tidal Force. The objective oF this 
thesis was to investigate these Findings by analyzing the 
gravitational Forces oF the moon and the sun with relation 
to recorded earthquakes. 

Two approaches in analyzing earthquake data are 
presented. The First approach was to compare the occurrence 
oF earthquakes with the Synodic Lunar Period, the Draconic 
Lunar Period, and the Anomalistic Lunar Period. The second 
approach was to investigate the tidal Force exerted by the 
sun and the moon with relation to 155 recorded earthquakes 
which occurred in the San Fernando Valley From April, 1971, 
to November, 1971. 
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The three lunar periods discussed in the first approach 
are as follows: 

1. Synodic Lunar Period, mean value 29.530589 solar 
days. Defined as new moon to new moon, where the 
new moon occurs when the geocentric longitudes of 
the sun and the moon are the same [Ref. 12] . 

2. Draconic Lunar Period, mean value 27.212220 solar 
days. Defined as the ascending node to the 
ascending node, ^^?here the ascending node pas- 
sage of the moon is from minus degrees latitude 
to positive degrees latitude [Ref. 12]. 

3. Anomalistic Lunar Period, mean value 27.554531 
solar days. Defined as perigee to perigee, 
where perigee is that point of the moon's orbit 

at which the moon is nearest to the earth [Ref. 12] . 

Five lunar periods [Ref. 12] actually exist but the 
above three were included for the reasons given below. 

In analyzing the earthquake data with respect to the 
Synodic Lunar Period, the occurrence of earthquakes with 
respect to the positions of the sun and the moon are 
examined. In analyzing the data with respect to the 
Draconic Lunar Period, the occurrence of earthquakes with 
respect to the moon's position in the Northern Hemisphere 
and Southern Hemisphere were examined. This may have been 
important since the data used included only earthquakes 
which occurred in California. In analyzing the earthquake 
data with respect to the Anomalistic Lunar Period, the 
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difference in the gravitational force of the moon on the 
earth was examined as the moon proceeded from perigee to 
apogee to perigee. 

The earthquake data received from the National Center 
for Earthquake Research spanned three time periods and 
occurred in two areas of California. The first set of data 
included 155 earthquakes recorded from April 1, 1971, to 
November 17, 1971. These are classified as aftershocks of 
the February 9, 1971, San Fernando Valley Earthquake 
(Figure 1). The magnitudes ranged from 1.0 to 4.2 as 
measured on the Richter Scale. This set of data will be 
referred to as the San Fernando Valley Data (SFVD) . 

The second set of data included 2879 earthquakes 
recorded from January 1, 1969, to August 23, 1970. These 
earthquakes occurred in the San Francisco Bay Area (Figure 
2) in the vicinity of the San Andreas Fault, the Hayward 
Fault, and the Calaveras Fault. The magnitudes ranged 
from 0.11 to 5.06 as measured on the Richter Scale [Refs. 

4 and 5] . This set of data will be referred to as the San 
Francisco Bay Area Data I (SFBAD-I). 

The third set of data included 2734 earthquakes recorded 
from August 1, 1970, to December 31, 1971. The area was 
the same as SFBAD-I. The magnitudes ranged from 0.29 to 
4.45 as measured on the Richter Scale [Refs. 5 and 6]. 

This set of data will be referred to as the San Francisco 
Bay Area Data II (SFBAD-I I). 
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Region (Indicated) of SFVD 



FIGURE 1 
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Region of SFBAD-I and SFBAD-II 
FIGURE 2 
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SFBAD-I and SFBAD-II both included earthquakes which 
occurred in August, 1970. In the analysis, as discussed 
in the next section, the duplicate earthquakes were eliminated. 

In analyzing the data with respect to the lunar periods, 
the following hypotheses were investigated: 

Null Hypothesis: Earthquakes occur in a 

uniform distribution with 
respect to the' lunar period. 
Alternative Hypothesis: Earthquakes do not 

occur in a uniform distri- 
bution vi^ith respect to the 
lunar period. 

Two tests were used in the analysis: The Chi-Square 

Test for Significance [Ref. 17] , and the Goodness of Fit 
Statistic Vn [Ref. 13] . The tests are discussed in v>:hat 
follows : 

The Chi-Square Test for Significance 

N, the total number of observations in the 
samples w'hich were analyzed. 

H, the total number of equal intervals in 
which the lunar period was divided. 

X = H , the expected number of observations v;hich 
occurred in each interval of a uniform 
distribution . 

Ti, the total number of observations recorded within 
each interval, i=l,...,H. 
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d£ = H - 1, the degrees of freedom associated with 
the test. 

H (T.-X)^ 

X , = — , the Chi-Square statistic. 

df j = i X 

The Goodness of Fit Statistic Vn: 

N, the total number of observations in the 
sample analyzed. 

FCx), the uniform distribution on the lunar 
period in which the origin was the origin 
of the lunar period. 

Fj^(x), the distribution of the observations on 
the lunar period with the same origin as F(x). 

X, the time of each earthquake in the sample. 

Vn = Sup (Fj^(x)-F(x)) - Inf (Fj^ (x) - F (x) ) , 

-oo<x<oo -oo<x<oo 

, the 

Vn statistic. 

Stephens [Ref. 13] stated that if the observations 
are points on a circle, the value of Vn obtained from the 
above equation does not depend on the choice of the origin 
for measuring x. The distribution of Vn refers to the dis- 
tribution under the null hypothesis [Ref. 13] . Tables of 
significance of the statistic /H Vn are included in 
Reference 13. 

A simulation of 1000 earthquakes over a period of six 
months was included in the analysis. The simulation was 
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performed to test the uniform distribution of the events 
over the three lunar periods. 

The final analysis in this thesis v\:as made possible 
by the use of a computer program by Professor Rex H. 

Shudde of the Naval Postgraduate School. This computer 
program calculated three components of the tidal force 
(measured in microgals) at the epicenter of an earthquake. 

The three components are the radial component, the 
north component, and the east component. The radial com- 
ponent of the tidal force vector is through the center of 
the earth and the epicenter of the earthquake. The north 
and east components of the tidal force vector are in the 
plane which passes through the epicenter and is perpendicular 
to the radial component. The computer program also calcu- 
lated the rates of change of the three components described 
above. The three components were plotted for a 24-hour 
time period on the day each earthquake occurred. A visual 
comparison of the plots of the 155 earthquakes which 
occurred in the San Fernando Valley with respect to the 
three components of the tidal force was made. The results 
of the comparison are in Section III. 



12 



II. ANALYSIS AND DISCUSSION OF THE LUNAR 



PERIOD COMPUTER PROGRAM 

The computer program used to analyze the frequency of 
earthquake occurrences during the Synodic, Draconic and 
Anomalistic Lunar Periods contains the following steps. 

The first step involves the time frame of the three sets 
of data and the lengths of the lunar periods within these 
time frames. The mean times of the three lunar periods, 
as given in the Introduction, were not used in the program, 
since the actual lengths of the periods vary considerably. 
An example of this difference was the Anomalistic Month 
which has a mean period of 27.32166 days. During 1971 the 
length of the Anomalistic Month ranged from 25.41667 
days to 28.45833 days [Ref. 9]. 

The times of new moon and lunar perigee for the years 
1969, 1970 and 1971 are tabulated in References 7, 8 and 
9. From these times (to be referred to as origins) the 
different period lengths of the Synodic and Anomalistic 
onths were computed. References 7, 8 and 9 also include 
tables indicating the apparent latitude of the moon for 
the respective years at intervals of one half day. The 
origins for the Draconic Month, as the moon passed from the 
Southern Hemisphere to the Northern Hemisphere, were 
computed from these tables by linear interpolation. This 
procedure should be fairly accurate since it involves only 
minutes and seconds of latitude. 
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The origins of the three lunar periods and the times of 
each earthquake occurrence were put into the program in 
year, month, day, hours, minutes and seconds in Greenwich 
Mean Time (GMT), Within the program the year, month and 
day of each event were converted to the respective Julian 
Day Number (JD) [Ref. 14], The hours, minutes and seconds 
constituted a fraction (Tract) of a day. The Julian Date 
was computed as follows: 

Julian Date = JD + Tract - 0.50 

The one half day was subtracted due to the differences 
in origins of the Julian Day Number and GMT. The Julian 
Day Number had its origin at noon whereas GMT has its origin 
at midnight. 

In order to simplify the statistical tests which followed 
in the program, the data was truncated to insure that no 
fractions of a period were included. An example of this 
was the STEAD- 1 and the Synodic iMonth. The time the 
STEAD- I started was January 1, 1969, which was between the 
two new moons of December 19, 1968, and January 18, 1969. 
Since the events that may have occurred from December 19, 
1968, to December 31, 1968, were not included, the events 
from January 1, 1969, to January 18, 1969, were excluded 
in the analysis. Similarly fractions of periods were 
eliminated for the Draconic and Anomalistic Months. 

Each of the three lunar time periods was reduced to the 
unit circle as explained in the following example. The 
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length of the Synodic Lunar Period from April 25, 1971, 
to May 24, 1971, was 29.354170 days. At the time of the new 
moon on April 25, the value on the unit circle was zero. 

At the time of the first quarter moon on May 2, the value 
on the unit circle was 0.25. At the time of the full moon 
on May 9, the value on the unit circle was 0.50. At the 
time of the third quarter moon on May 17, the value on the 
unit circle was 0.75. Finally at the time of the new moon 
on May 24, the value on the unit circle was 1,0 [Ref. 9], 

The time at which each earthquake occurred during the 
lunar period (residual) was placed upon the unit circle 
with respect to the respective origins. This was accom- 
plished by dividing the difference between the Julian Date 
of the earthquake and the Julian Date of the origin immedi- 
ately preceding the earthquake by the length of the period 
in which the earthquake occurred. 

The unit circle was divided into 20 equal intervals 
in order to perform the Chi-Square Test for Significance 
[Ref. 17] as explained in the Introduction. Ostle [Ref. 10] 
suggests, as a rule of thumb, to use no fewer than three for 
the expected frequency in each interval. Ostle [Ref. 10] 
further contends that if some expected numbers are too small 
the Chi-Square statistic will be a poor indicator of the 
validity of the hypothesis under test. By selecting 20 
intervals, the earthquakes in the three sets of data comply 
with the rule of thumb for expected frequencies. 
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The residuals that were computed were then sorted into 
ascending order [Ref. 2] so that the Vn Goodness of Fit 
Test could be computed as explained in the Introduction. 

The Vn Goodness of Fit Test is a modification of the 
Kolmogorov- Smirnov Test [Ref. 13] which is considered to 
be generally a more powerful test than the Chi-Square 
Test [Ref. 10]. 

The output of the program consisted of the total number 
of events that occurred during each of the three lunar 
periods; the number of events and the percentages whicli 
occurred within each of the 20 intervals, the Chi-Square 
statistic, and the Vn and /M Vn statistics. 
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III. RESULTS 



The simulation of the 1000 earthquakes as stated in the 
Introduction was performed to test the uniform distribution 
of the events on the three lunar periods. The null hypothe- 
sis that the simulated earthquakes were distributed uniform- 
ly was accepted for each of the lunar periods. The bases 
for the acceptances were the Chi-Square Test with 19 
degrees of freedom and the Vn Test. The Chi-Square statis- 
tics and the Vn statistics calculated are given below: 





Synodic 


Draconic 


Anomalistic 


. Chi-Square 


23.28 


22.04 


28.96 


/N Vn 


1.4002 


1.2130 


1.3995 



The null hypothesis and alternative hypothesis are: 

Null Hypothesis: Earthquakes occur in a 

uniform distribution with respect 
to the lunar period. 

Alternative Hypothesis: Earthquakes do not occur 

in a uniform distribution with 
respect to the lunar period. 

If the null hypothesis is accepted, it might imply that 
the moon's gravitational effects on the earth do not trigger 
earthquakes. Conversely if the null hypothesis is rejected 
and the alternative hypothesis is accepted then it might 
imply that the moon's gravitational effects on the earth 
do trigger earthquakes. 
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The results of the statistical tests performed on the 
three sets of data with respect to the Synodic, Draconic 
and Anomalistic Lunar Periods are presented in tables. 
Results suggested by the tidal force plots are discussed 
following the tables. 

In each table the following definitions are used: 

Mag - Magnitudes of earthquakes greater than or 

equal to the magnitude entered in the table. 

No - Number of earthquakes in the sample. 

Chi - The Chi-Square statistic with 19 degrees 
of freedom. 

a - The one per cent significance level for the 
Chi-Square statistic [Ref. 11]. 

Vn - The Vn statistic. 

Vn/N - The Vn statistic for large values of No. 

b - The one per cent significance level for the Vn 
statistic . 

Acc - Acceptance of the null hypothesis at the one 
per cent significance level. 

Rej - Rejection of the null hypothesis at the one 
per cent significance level. 

Polar plots of the occurrence of earthquakes during the 
lunar periods are presented in Figures 3, 4, 5, 6, 7, 8, 9, 
10, and 11. 
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A. SAN FERNANDO VALLEY 



1. Synodic Month 





Mag 


No 


Chi 


a 


Vn 


Vn/N 


b 




All 


110 


43.818 


REJ 


.18258 


1.91496 


ACC 


> 


1.75 


106 


33.245 


ACC 


.15788 


1.62552 


ACC 


> 


2.00 


104 


34.077 


ACC 


.15444 


1.57496 


ACC 


> 


2.25 


96 


38.167 


REJ 


.14963 


1.46607 


ACC 


> 


2.50 


88 


28.364 


ACC 


. 14395 


1.35035 


ACC 








Table 


1 








2. 


Draconic Month 












Mag 


No 


Chi 


a 


Vn 


Vn/N 


b 


> 


All 


119 


34.109 


ACC 


.16705 


1.82225 


ACC 


> 


1.75 


115 


26.739 


ACC 


.14775 


1.58449 


ACC 


> 


2.00 


112 


25.857 


ACC 


.14285 


1.51158 


ACC 


> 


2.25 


103 


28.068 


ACC 


.14267 


1.44793 


ACC 


> 


2.50 


93 


26.355 


ACC 


.14475 


1.39593 


ACC 



Table 2 
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3. Anomalistic Month 





Mag 


No 


Chi 


a 


Vn 


Vn/N 


b 




All 


122 


21.934 


ACC 


.15824 


1.74781 


ACC 


> 


1.75 


118 


19.627 


ACC 


.14189 


1.54135 


ACC 


> 


2.00 


115 


20.478 


ACC 


. 13843 


1.48455 


ACC 


> 


2.25 


104 


21.000 


ACC 


.14804 


1.50976 


ACC 


> 


2.50 


93 


19.903 


ACC 


.14287 


1.37782 


ACC 



Table 3 
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FIGURE 5 
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SAN 


FRANCISCO BAY 


EARTHQUAKE 


: DATA 


I 






1. 


Synodic 


Month 














Mag 


No 


Chi 


a 


Vn 


Vn/R 


b 




All 


2721 


156.405 


REJ 


.05351 


2.79130 


REJ 


> 


0.50 


2554 


171.090 


REJ 


. 05848 


2.95563 


REJ 


> 


0.75 


2333 


187.660 


REJ 


.06551 


3.16441 


REJ 


> 


1.00 


2008 


202.358 


REJ 


.06693 


2.99933 


REJ 


> 


1.25 


1580 


204 . 278 


REJ 


.07694 


3.03841 


REJ 


> 


1.50 


1193' 


186.966 


REJ 


.08954 


3.09260 


REJ 


> 


1.75 


803 


170.425 


REJ 


.09327 


2.64308 


REJ 


> 


2.00 


483 


106.772 


REJ 


.09356 


2.05609 


REJ 


> 


2.25 


301 


93.485 


REJ 


.11281 


1.95726 


ACC 


> 


2.50 


177 


125.034 


REJ 


.16622 


2.21147 


REJ 


> 


2.75 


111 


98.189 


REJ 


.18296 


1.92756 


ACC 


> 


3.00 


64 


83.500 


REJ 


.24541 


1.96325 


REJ 



Table 4 
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2721 EARTHQUAKES 
PERCENTAGES INDICATED 



FIGURE 6 



2. Draconic Month 



Mag 

All 

> 0.50 

> 0.75 

> 1.00 

> 1.25 

> 1.50 

> 1.75 

> 2.00 

> 2.25 

> 2.50 

> 2.75 

> 3.00 



No 


Chi 


a 


Vn 


Vn/N 


b 


2780 


127.611 


REJ 


. 04939 


2.60438 


REJ 


2611 


141.823 


REJ 


.04908 


2.50802 


REJ 


2384 


162.342 


REJ 


.05230 


2.55383 


REJ 


2 04 8 


163.230 


REJ 


.06090 


2.75592 


REJ 


1611 


180.521 


REJ 


.06805 


2.73115 


REJ 


1214 


148.965 


REJ 


.07410 


2.58155 


REJ 


818 


115.203 


REJ 


.08995 


2.57271 


REJ 


493 


65.823 


REJ 


. 09078 


2.01568 


REJ 


310 


48.193 


REJ 


.11873 


2.09050 


REJ 


184 


70.348 


REJ 


.14956 


2.02866 


REJ 


116 


58.138 


REJ 


.13858 


1.49258 


ACC 


66 


45.515 


REJ 


.14618 


1.18757 


ACC 



Table 5 
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NODE 
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FIGURE 7 
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3. Anomalistic Month 



Mag 

All 

> 0.50 

> 0.75 

> 1.00 

> 1.25 

> 1.50 

> 1.75 

> 2.00 

^ 2.25 

> 2.50 

> 2.75 

> 3.00 



No 


Chi 


a 


Vn 


Vn>/N 


b 


2794 


148.190 


REJ 


.05336 


2.82060 


REJ 


2625 


144.379 


REJ 


.05320 


2.72548 


REJ 


2398 


153.409 


REJ 


. 05437 


2.66237 


REJ 


2061 


156.613 


REJ 


.06303 


2.86128 


REJ 


1623 


142.878 


REJ 


.06268 


2.52532 


REJ 


1224 


115.314 


REJ 


.07092 


2.48117 


REJ 


822 


94.448 


REJ 


.08887 


2.54783 


REJ 


496 


51.823 


REJ 


.09322 


2.07606 


REJ 


310 


40.710 


REJ 


.10603 


1.86686 


ACC 


184 


38.826 


REJ 


.12593 


1.70822 


ACC 


116 


34.345 


ACC 


.17119 


1.84375 


ACC 


66 


26.727 


ACC 


.14634 


1.18887 


ACC 



Table 6 
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APOGEE 




SFBAD-I 

ANOMALISTIC MONTH 
2794 EARTHOUAKES 
PERCENTAGES INDICATED 



FIGURE 8 



SAN 


FRANCISCO BAY 


EARTHQUAKE DATA 


II 






1. 


Synodic 


Month 














Mag 


No 


Chi 


a 


Vn 


Vn/N 


b 




All 


2595 


135.767 


REJ 


.06405 


3.26256 


REJ 


> 


0.50 


2571 


132.034 


REJ 


.06204 


3.14580 


REJ 


> 


0.75 


2438 


132.336 


REJ 


.06392 


3.15630 


REJ 


> 


1.00 


2144 


122.847 


REJ 


.06422 


2.97360 


REJ 


> 


1.25 


1691 


96.002 


REJ 


. 06416 


2.63846 


REJ 


> 


1.50 


12 2 8 


78.384 


REJ 


.07451 


2.61101 


REJ 


> 


1.75 


810 


65.506 


REJ 


.08914 


2 . 53705 


REJ 


> 


2.00 


459 


47.710 


REJ 


.11265 


2.41349 


REJ 


> 


2.25 


275 


35.764 


ACC 


.11716 


1.94280 


ACC 


> 


2.50 


155 


33.774 


ACC 


.14005 


1.74338 


ACC 


> 


2.75 


79 


27.076 


ACC 


.18849 


1.67538 


ACC 



Table 7 
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FULL MOON 




NEW MOON 



FIGURE 9 



2. Draconic Month 



Mag 

All 

> 0.50 

> 0.75 

> 1.00 

> 1.25 

> 1.50 

> 1.75 

> 2.00 

> 2.25 

> 2.50 

> 2.75 



No 


Chi 


a 


Vn 


Vn/F? 


b 


2556 


31.512 


ACC 


.03189 


1.61201 


ACC 


2534 


32.140 


ACC 


.03155 


1.58830 


ACC 


2406 


31.506 


ACC 


.03599 


1.76540 


ACC 


2117 


23.680 


ACC 


.03257 


1.49860 


ACC 


1668 


17.875 


ACC 


.03807 


1.55467 


ACC 


1212 


25.294 


ACC 


.05430 


1.89038 


ACC 


810 


18.151 


ACC 


.05296 


1.49892 


ACC 


454 


26.529 


ACC 


.06228 


1.32712 


ACC 


271 


26.491 


ACC 


. 07411 


1.22006 


ACC 


152 


31.158 


ACC 


.09405 


1.15959 


ACC 


77 


25.597 


ACC 


.13345 


1.17102 


ACC 



Table 8 
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SFBAD-II 
DRACONIC MONTH 
2556 EARTHQUAKES 
PERCENTAGES INDICATED 



FIGURE 10 



b 



3. Anomalistic Month 





Mag 


No 


Chi 


a 


Vn 


Vn/N 




All 


2634 


36.752 


REJ 


.04930 


2.53041 


> 


0.50 


2611 


35.181 


ACC 


.04856 


2.48153 


> 


0.75 


2479 


40.943 


REJ 


.05378 


2.67761 


> 


1.00 


2184 


45.982 


REJ 


.05309 


2.48120 


> 


1.25 


1717 


43.641 


REJ 


.05110 


2.11722 


> 


1.50 


1251 


48.265 


REJ 


.06704 


2.37127 


> 


1.75 


829 


49.673 


REJ 


.07730 


2.22563 


> 


2.00 


473 


36.302 


REJ 


.09228 


2.00702 


> 


2.25 


286 


31.203 


ACC 


. 10886 


1.84091 


> 


2.50 


162 


32.321 


ACC 


.11470 


1.45990 


> 


2.75 


84 


31.238 


ACC 


.14562 


1.33464 
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SFBAD-II 

ANOMALISTIC MONTH 
2634 EARTHOUAKES 
PERCENTAGES INDICATED 



FIGURE II 
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D. TIDAL FORCE COMPONENTS WITH RESPECT TO THE 

SAN FERNANDO VALLEY DATA 

Figure 12 is an illustration of the tidal force plots 
used in the analysis. The occurrence of the earthquake 
in the plot was May 19, 1971 at 0451 GMT. The magnitude 
of the earthquake was 2.6. The radial component of the 
tidal force, as measured in microgals, was positive and 
decreasing in value (negative slope). The north component 
was negative and increasing in value (positive slope) . 

The east component was negative and decreasing in value 
(negative slope) . 

The same 110 SFVD earthquakes which were analyzed with 
respect to the Synodic Lunar Period were investigated as 
the event presented above. The results are tabulated in 
Table 10. The time frame in which these earthquakes 
occurred was from April 25, 1971, to October 19, 1971. 

Ten thousand earthquakes uniformly distributed in time 
were simulated during the previously mentioned time frame 
and in the same area as the SFVD. These simulated earth- 
quakes were analyzed in the same manner. The results are 
also tabulated in Table 10. 

Each component at the time of the earthquake was either 
positive or negative in value. The derivatives of each 
of the three components was either positive or negative 
at the time of the earthquake. Sixty-four combinations of 
the three components and their slopes were then considered. 
In Table 10 the following definitions are used. 
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figure 12 



RAD - 


The radial component of the tidal force. 


NOR - 


The north component of the tidal force. 


EAS - 


The east component of the tidal force. 


REC - 


The number of earthquakes recorded with the 
same combination. 


PREC - 


The percentage of the earthquakes recorded 
with the same combination. 


SIM - 


The number of simulated earthquakes recorded 
with the same combination. 


PSIM - 


The percentage of the simulated earthquakes 
with the same combination. 


+ 


Component of tidal force was positive. 
Component of tidal force was negative. 


P 


Value of the component was increasing. 


n 


Value of the component was decreasing. 



The data in Table 10 were used to compare the distribu- 
tion o£ the earthquakes with the distribution of the simulated 
earthquakes by means of a Contingency Table [Ref. 10], 

The null hypothesis was that the two distributions are the 

\ 

same. The combinations of the three components of the 
tidal force vector and the three derivatives, for which 
the actual and the simulated earthquakes never occurred, 
were eliminated. This resulted in a 2X34 Contingency Table. 
The Chi-Square statistic with 53 degrees of freedom was 
computed (39.972) and the null hypothesis was accepted. 
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REC 


PREC 


SIM 
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4.55 
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21 
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+P 


4 
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4.69 


+n 
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+P 
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-n 


-P 


+P 


0 


0 


1 


.01 


+P 


+n 


+P 


0 


0 


74 


. 74 


-P 


+n 


+P 


3 


2.73 
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1,79 


+n 


+n 


+P 


0 


0 


0 


0 
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+n 


+P 


0 


0 


0 


0 


+P 


-n 


+P 


1 


.91 


215 


2.15 


-P 
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+P 
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8.18 
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0 


0 


0 
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69 
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-P 
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4.55 
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0 


0 
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-P 


0 


0 


17 


.17 


+n 


-P 


-P 


4 


3.63 
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2.46 


-n 


-P 


-P 


14 


12.73 


1356 


13.56 


+P 


+n 


-P 


0 


0 


0 


0 
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RAD 


NOR 


HAS 


REC 


PREC 


SIM 


PSIM 


-P 


+n 


"P 


0 


0 


0 


0 


+n 


+n 


-P 


0 


0 


0 


0 


-n 


+n 


-P 


1 


.91 


236 


2.36 


+P 


-n 


'P 


0 


0 


0 


0 


-P 


-n 


’P 


1 


.91 


5 


.05 


+n 


-n 


"P 


0 


0 


23 


.23 


-n 


-n 


-P 


6 


5.45 


479 
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+P 
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*^P 


+n 
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+n 
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+n 


0 


0 
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0 
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1.06 
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+n 
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3.63 
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0 


0 


21 
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0 


0 


5 
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0 
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0 


19 
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+n 


-n 


+P 


-n 


-n 


-P 


-n 


-n 


+n 


-n 


-n 


-n 


-n 


-n 



PREC 


SIM 


PSIM 


0 


0 


0 


0 


13 


.13 


0 


14 


.14 


7.27 


1617 


16.17 


7.27 


579 


5.79 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


.01 


0 


0 


0 


.91 


97 


.97 


.91 


69 


.69 



REC 

0 

0 

0 

19 

8 

0 

0 

0 

0 

0 

0 

1 

1 
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IV. DISCUSSION OF RESULTS AND CONCLUSIONS 



The occurrence of earthquakes, at the larger sample 
sizes, was rejected as a uniform distribution with respect 
to the lunar periods, in all cases except one (Table 8). 

On certain days within the time frames of SFBAD-I and 
SFBAD-II, a larger number of earthquakes occurred relative 
to other days in these time periods. This "sv^:arming'' 
effect of numerous earthquakes occurring within hours and 
minutes of each other may indicate the lack of independence 
between the events. The "swarming” effect could also indi- 
cate that earthquakes do occur in a non-uniform manner. In 
order to decrease the dependence between the events, if 
the dependence does exist, the data were analyzed according 
to the magnitudes of the earthquakes. This method elim- 
inated the "swarming" effect. An example of the "s^^farming" 
effect was on June 10, 1969, v\?hen 15 earthquakes occurred. 
When magnitudes greater than or equal to 2.0 were examined, 
the number of earthquakes which occurred on June 10, 1969, 
was reduced to five. The null hypothesis was still rejected 
for the three lunar periods. The fact remains that the 
null hypothesis was rejected in all cases but one at the 
larger sample sizes. 

It is interesting to compare the polar plots of the 
Synodic Lunar Periods of the SFVD and the SFBAD-II. The 
largest percentage of earthquakes occurred, in both cases, 
in the first interval, or at the time of the new moon. 
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Also, the time periods of the two sets of data were similar. 
However, the polar plot of the SFBAD-I was not similar, but 
neither was the time period. 

The occurrence of earthquakes with respect to the 
Draconic Lunar Period in the SFBAD-I and SFBAD-I I are not 
the same. When analyzing the SFBAD-I, the null hypothesis 
was rejected at all magnitudes, but in analyzing the 
SFBAD-II the null hypothesis was accepted at all magnitudes. 

The occurrence of earthquakes with respect to the Anoma- 
listic Month in the SFBAD-I and SFBAD-II are similar to 
each other. The rejection of the null hypothesis at the 
one per cent significance level was the same for the larger 
sample sizes. 

The null hypothesis was rejected when analyzing the 
SFVD with respect to the Synodic Lunar Period at the one 
per cent significance level when using the Chi-Square Test 
but was accepted using the Vn Test (Both tests would have 
rejected the null hypothesis at the 2.5 per cent signifi- 
cance level). The rejection was due mainly to the occur- 
rence of 17 of the 110 earthquakes near the time of the 
new moon. In comparing these 17 earthquakes with respect 
to the combinations of the three components of the tidal 
force (as discussed in Section III), nothing of interest 
seems to be revealed. Four of the 17 earthquakes had the 
same combination but the four also occurred on the same 
day. The remaining 13 earthquakes had different combinations. 
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The acceptance in Section III, that the distributions 
of the 110 actual earthquakes and the 10,000 simulated 
earthquakes were the same seemed quite significant. Further 
analysis into why earthquakes never occur with certain 
combinations and why large numbers of earthquakes occur with 
other combinations may prove or disprove that the tidal 
force can trigger earthquakes. 
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V. THE COMPUTER PROGRAM 



The time (GMT) of each origin of the lunar periods 
was read into the computer by the year (NA) , the month (NB) , 
the day (NC) , the hour (ND) , the minute (NE) , and the 
second (SE) . As discussed in Section II, the time of each 
origin was converted [Ref. 14] to the corresponding Julian 
Date (ORGP(I,J)). 

The time (GMT) of each earthquake was read into the 
computer by the year (MA) , the month (MB) , the day (MC) , 
the hour (MD) , the minute (ME) , and the second (SE) . As 
discussed in Section II, the time of each earthquake was 
converted [Ref. 14] to the corresponding Julian Date 
(JDGMT(I)). Subsequently each of the JDGMT(I) v^^as con- 
verted to the residual (RESD(I,J)) for each lunar period. 

The following list defines the terms used in the Lunar 
Period Computer Program; 

NBOX The number of equal intervals on the 



JDAY 



0RGP(1,I) 

0RGP(2,I) 

0RGP(3,I) 

PER(1,I) 

PER(2 ,1) 

PER(3,I) 



interval (0,1) ^ 

t h 

The origin of the I Synodic Month. 

t ll 

The origin of the I^ Draconic Month. 

The origin of the I^ Anomalistic Month. 

The period of the I^^^ Synodic Month. 

*t ll 

The period of the I^ Draconic Month. 

ll 

The period of the I^ Anomalistic Month. 
Subroutine to convert year, month and day 
to the Julian Day Number [Ref. 14] . 
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BEG(l) 


The first origin of the Synodic Month 
to occur within the time frame of the data. 


BEG (2) 


The first origin of the Draconic Month to 
occur within the time frame of the data. 


BEG (3) 


The first origin of the Anomalistic Month 
to occur within the time frame of the data. 


EN(1) 


The last origin of the Synodic Month to 
occur within the time frame of the data. 


EN(2) 


The last origin of the Draconic Month to 
occur within the time frame of the data. 


EN(3) 


The last origin of the Anomalistic Month 
to occur within the time frame of the data. 


JDGMT 


The Julian Date each earthquake occurred. 


JDLP 


The time which each earthquake occurred 
within the period. 


RES 


The residual computed for each earthquake 
and placed on the unit circle. 


NCNT 


The number of occurrences within each 
sub interval . 


NPTS 


The total number of earthquakes. 


EXVAL 


The expected number of earthquakes within 
each subinterval. 


PCT 


The percentage of the number of occurrences 
within each subinterval. 


CHISQ 


The Chi-Square statistic. 


SHELL 


Subroutine to sort the residuals in 
ascending order [Ref. 2]. 
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DIFF 



The Vn statistic. 



ROOTS 



The Vn statistic times /H. 
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onnoo ooooo 



0001 



7000 

^ V ^ 



*1 ^ ^ 



6999 



6005 

6000 



6995 

6015 



0040 

0030 



0045 

0035 

0005 

1500 

%JU «jU ^ 
^ n' -r V 






0020 



REALMS FRACT, JDGMT, JDLP , RES , PER, PE RP.ORG ,0RGP 
DIMENSION PERI 5 ) ,GRG(5 ) , JDLP (5 ) , NCNT ( 30 , 5 ) , RESD { 1 60 , 5 ) 
DIMENSION CHISQ(5) , SUM { 5 ) , EM I N { 5 ) , EM AX ( 5 ) , D I F F { 5 ) 
DIMENSION ROOTS! 5) ,ORD( 161 ) 

DIMENSION ORGP! 3 , 1 0 ) , P ERP ( 3, 9 ) , PCT { 30, 5 ) , EXVAL (5) 
DIMENSION NUM(5) ,JDGMT(160) ,3EG(3) ,EN{3), ENPTS(5) 

DATA CHlSQ/5=f^0.0/, SUM/5^-0. 0/ 

DATA NUM/5=^0/ 

REA0(5,0001) NBOX 
FORMAT (14) 

BOX=NBGX 

WRITE! 6 ,7000) NBOX 

FORMAT !• 1 ',/, 42X, 'UNI FORM ?',I5,‘ BOXES',/) 

CON\/ERT ORIGIN GREGORIAN DATES TO JULIAN DATES 

OU OU JU JU JU ^ ^ ^ 

-"j-* ^ ^ ^ ^ ^ Y ^ ^ ^ ^ ^ V *1^ ^ ^4^ 'ip 'P ^ ^ ^P ’Hr 'P *P ^P ^ 

DO 6000 1=1,3 
DO 6005 J=1 , 10 

READ! 5,6999) NA , NB , NC , ND , NE , SE 
FORMAT! IX, 14,13, 13, 14, 13, F6. 2) 

CALL J0AY!NC,N3,NA, JD) 

ORGP ! I , J )= JD+! FLDAT(ND)+FLOAT! NE) /60 . 0+S E / 3600. 0 ) /24 . 0 

ORGP ! I , J ) = ORGP ( I , J )-0. 50 

CONTINUE 

CONT INUE 

BEG! 1 )=ORGP! 1, 2 ) 

BEG!2) =ORGP(2,2) 

BEG! 3)=0RGP! 3, 2) 

EN(1)=CRGP!1,3 ) 

EN! 2) =0RGP! 2,9 ) 

EN!3 )=0RGP( 3, 9 ) 

ORG(4) =ORGP!2, 1 ) 

ORG! 5) =ORGP! 2, 1 ) 

DO 6015 1=1,3 

READ(5,6995) ( P ERP ! I , J ) , J = 1 , 9 ) 

FORMAT! 10F8.5) 

CONT INUE 

PER! 4) =27. 32166 

PER(5)=27.32158 

NPTS=0 

DO 0030 1=1, NBOX 
DO 0040 J= 1, 5 
NCNT!I , J)=0 
CONT INUE 
CONT INJE 
DO 0035 1=1,160 
DO 0045 J=l,5 
RESD! I , J ) = 0.0 
CONTINUE 
CONT INUE 
11=1 

RE AD ! 5 , 1 5 00 , EN D = 002 0 ) MA , MB , MC , MD , ME , S E 
FORMAT! IX, 14, I 3, 13, 14, I 3, '=6. 2) 

IF(MA. EQ.O ) GO TO 0020 

U- -u -JU -i- ^ .JO U. V- ■*' ^ V Jr J- 'R' ■J' J- V- *r^ J' 'J' -V ^ ^ — -J- ^ 

T- «r* *r* ip ... ^ •r' T' -T* •V *r* 'T- TT* T' *1- V* 'T' 'T' ^ pp pp .p .p. .p pp r,. pp ^ ^r- nr 



CONVERT EARTHQUAKE GREGORIAN DATES TO JULIAN OATES 



OU OU 4ju OU wV ^ w 4.1. «JU ou ^ ^ Ju JU sV 4/4 JU JU *X4 w OU JY ^ ^ 

nr^ ^ ^ ^ ^ ^ ^ Y. ^ ^ ..^N ^ ^ ^ ^ ^ ^ :jc ... ^ ^ ^ ^ # h ^ 

CALL JDAY! MC.MB, MA, JD) 

NPTS=NPTS+1 

FRACT= ( FLOAT! MD)+FLOAT (ME) /60. O+SE/3600.0) /24.0 

JDGMT! NPTS ) = JD+F RACT-0 . 500 

GO TO 0005 

DO 4000 1=1,5 

DO 4005 L=1,NPTS 

IF!I.GT.3) GO TO 4100 

IF! JDGMTIL ) .LI .BEG! I ) ) GO TO 4005 

IF! JDGMT(L) .GT.EN! I) ) GO TO 4005 

DO 6030 J=2, 10 
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I F( JDGMT (L ) .Gr .ORGP( I , J ) ) GO TO 6030 
ORG(I) =ORGP( It J-1 ) 

PER( I )=PERP( I , J-1) 

J=10 

6030 CONTINUE 

4100 NUM( I ) = MUM (1)4-1 

JDLP( I ) =JDGMT ( L)-ORG( I ) 

C *JL. ^ '<kV ^ tJU tJ> ^ ^ ^ »'-• ^ ^ ^ ^ V* ^ ^ ^ JU w Vt» ^ ^ -Jf ^ JU aJU w ^ <A« JU 

'r* -r- ' *T -Y- ¥■ -T- <v v ' — r ■4' 'i' h' ^ - 4' "r Y* 'r ^ ^ ,r- ^ 'r ^ ^ *t' ^ ^ ^ ^ ^ ^ ^ ^ .*r ^ ^ ^ ^ ^ ^ ^ ^ ^ 

C 

C COMPUTE RESIDUALS 

C 

R5S=Dt10D( JDLP( I),PER(I))/PER(I) 

IF(RES. LT. O.DO) RE S =RE S+- 1 . DO 

^ ^ ipp ^ >P ^ 4p ^ ^P «tp ^ ^ .^1% >|S ^ ^ -TP ^ #!*• ^ jp *p #pt ^ *rp -►P ^p #P ^ ^ ip <ip -p 

C 

C PLACE RESIDUALS INTO 20 BOXES ON INTERVAL (0,1) 

C FOR EACH LUNAR PERIOD 

C 

Q >^ :J? ^c * <c ^ >Jc *r *:>; 3)c # ^ :{t *# 5^ jjt * s^t tjr *# -c # 

INDEX=PES--;=BOXf 1.0 
NCNTdNDEX, I)=NCNT( INDEX, D-H 
RESD( I I , I )=RES 
I I = I I-H 
4005 CONTINUE 

ENPTS( I )=NUM( I ) 

EXVAL( I ) = ENPTS ( I )/BOX 
DO 0070 J=1,N30X 

PCT( J, I ) = FL0A1 (NCNT( J, I ) ) /ENPTS( I ) 

0070 CONTINUE 

C s!cs!!c:?'*5lr:&*«iJc;X5iiric^iic5ic**;ic:Jc4c4c5kic3k*s!:jk:ii5kJktir:i::fe&:icrfc5*^5k*:risric^5k 

c 

C COMPUTE CHI-SQUARE ON BOXES 

C 

C *J^ %V V -» %V « ^ >!•*' «»rV Ot#* %V v<^ ^ ^ ^ ^ ^ ^ «V 

rfp >p ^p P’*» ^ % ^p ^ ■» >P pF^ -p p- ppH >p pp p^ pt» ^p pp /p 4»p >p «*»% #p ^p p\ pp #p 

DO 0098 J=1,NB0X 

SUM( I ) = ( FLOAT! NCNT ( J ,I ) ) -EXVAL ( I ) ) *=i=2/ EXV AL ( I ) 

CHISQ( I )=CHIS0 ( I )-^SU^U I ) 

0098 CONTINUE 

C «JU ^ ^ ^ ^ w V JU ^ Ji# JU ^ ^ ^ ^ JU ^ JU ^ ^ JU ^ ^ ^ ^ JU JU ^ %U ju JU JU u# ^ ^ ^ ^ 

^ pp ^p »p pp pp >p »p "p n*v p -p ijp up p% ^p #p jp ^|p ^p pp p% p pp p« pp #p p -p pp p« p-« >p up up pp pp pp up A% pN p-^ .p ^p 

c 

C SORT RESIDUALS 

C 

C J^ JU V-^ ^ «U JU JU J«* JU JU ^U JU Ji# J«« Ju JU JU %U JU JU JU JU JU JU JU JU JU J.« JU J-r JU J^ Jt» JU >U JU JU JU JU JU JU JU JU JU J-f JU JU JU JU JU JU 

pi^ ^ pr -'A ^ P-* P -r^ P '*** ^ ^^ ^'^ ***^ n -»♦ p* p ^ ^ nr nr ^ p% n**. p* ^ -^p ip n^ ^ ^ pp ^ ;p ^ ^ up ^ ^ jp pp n-^ n% ?P p% ^ p 

CALL SHELL! NUM ( I ) , RESD ( 1, I ) ) 

N0=NUM( I ) 

DO 0033 J=1,N0 

0RD( J)=FL0AT( J-1 )/ENPTS( I ) 

0033 CONTINUE 

0RD(NJM( I )+l)=1.0 
EMIN( I ) =0. 0 
EMAX( I ) =0. 0 

c 

C COMPUTE MODIFIED K-S TEST ON THE LUNAR PERIODS 

C 

C JU Ju JU JU JU J-» JU JU JU JU Ji# J^ JU JU A* JU JU JU J^ J** ^U JU JU JU JU JU Ju %V JU JU >V JU JU JU JU «Av JU JU JU «kA^ J^ %U JU %V JU «V J- JU 
n^ pp* *1^ p% p-» ♦r* -p V* nr *•“* nr nr* nr ^ V* nr n^* nr nr ^ ^ n** nr nr nr nr n* n^ nr nr n** nr n”^ n** ^ np ^ pp up u<* ,ij. up up .r-. pp .f**. 

DO 0054 J=1,N0 
RESDX=RES0 ( J , I ) 

EMAX( I ) =AMAX1( RESDX-0RD( J) ,EMAX( I ) ) 

EMIN ( I )=AM INK RE SD X-ORD ( J -H ) , EMIN( I ) ) 

0054 CONTINUE 

DIFF ( I ) = EMAX( I ) -EM INI I ) 

ROOTS ( I )=DIFF( I )MSQRT ( ENPTSI I ) ) 

11=1 

4000 CONTINUE 

C J«* ^4* JU J« <J^ %V ^ Jp J« JU J^- JU JU Jp JU JU JU JU JU ^U J-» JU JU JU JU JU JU JU JU JU JU JU w JU JU JU JU JU JU JU JU JU JU JU 

pp«t^ «p »p» pp ^ . pk rp n** ^ • n^ ^v* ^ ^ ^ np n** n’* n** n^ nr ^** n^ n* *r* n^ ^ n^ •r* nr nr n* n^ nr nr * ^ 

c 

C TABULATE THE RESULTS 

C 

'T'r ^ ^ ^ .f, 0^^ 0^ 0^ 0^ 0^ 0f^ 0^ 0^ 0,^ 0f^ ^ 0^ 0,^ 0^ 0^, ^ 0Y^ 'i' 'i' 

WRITE! 6,9999) 
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9999 FORMAT (9X, ' SYM3DIC MON TH ' , 5X , • DR ACOM I C MONTH', 5X, 

1 ' ANOMALI STIC MONTH ' , 5X , ' S I DER EAL MONTH', 5X, 

2'TROPICAL MONTH' ,//) 

DO 0060 I=1,NB0X 
WRITE(6,9900)( NCNT( I,J) ,J = 1,5) 

9900 FORMAT! 1 5X , I 4 , 1 X , I 4 , 1 3X , I 4 , 1 5X , I 4 , 1 5X , I 4 ) 

0060 CONTINUE 

WRITE (6,4999) ( NUM( J) ,J=1 ,5) 

4999 FORMAT ( / , 1 IX, I 4, ' QUAK E S ' , 8X , I 4 , ' QUAKE S ' , 1 1 X , I 4 , 

1' QUAKES' ,8X, I 4, ' QUAK 6S ' , 8X , 1 4 , ' QUAKES',/) 

WRITE! 6,7005) 

7005 FORMAT ! //,46X, 'PERCENTAGES' , /) 

WRITE (6,9999) 

DO 0155 I=1,N30X 
WRITE!6,8100)!PCT! I,J) ,J = 1,5) 

8100 FORMAT! 11X,F7.3 ,12X,F7.3 ,12X,F7.3,15X,F7.3, 12X, F7.3) 
0155 CONTINUE 

WRITEI6 ,2044 ) 

2044 FORMAT! ' 1' ,// ,45X, 'CHI -SQUARE TEST',/) 

WRITE! 6,9999) 

WRITE!6,3200)!CHISQ!J),J=1,5) 

8200 FORMAT! 11X,F7.3,12X,F7.3,12X,F7.3,15X,F7.3,12X, 

1F7. 3, //////////) 

WRITE! 6 ,9999) 

WRITE! 6,4988) ! NUM! J ) ,J=1 ,5) 

4988 FORMAT ! 2X, 9X, I 4, ' QUAK ES ' , 8X , I 4, ' QUAKE S ' , 1 1 X , I 4 , 

1' QUAKES' , 8X, I 4, ' QUAK ES ' , 8X , 14 , ' QUAKES',//) 
WRITE!6,2022) ( DIFf=! J) ,J = 1,5) 

2022 FORMAT ! 1 X , ’ MC D I F I E D ' , / , 3X , ' K-S ' , 3X , F 9. 5 , 1 0 X, F 9 . 5 , 1 OX , 
1F9.5,13X,F9.5,10X,F9.5,/,1X,'STATISTIC',//) 
WRITE!6,2033)( ROOTS! J) ,J=1,5) 

2033 FORMAT !2X, ' AdOV E ' , / , 1 X , 'STATISTIC',/,2X,'TIMES',2X, 
1F9.5,10X,F9.5,10X,F9.5,13X,F9.5,10X,F9.5,/,1X, 

2 ' ROOT ! NPTS ) ' ) 

STOP 

END 



subroutine JDAY! D,M,Y, j ) 

INTEGEF;*4 D,Y,YA,C 
IF!M.GT.2) GO TO 0055 
IM=M+9 
I Y=Y-1 
GO TO 0065 
0055 IM=M-3 
I Y=Y 

0065 C=IY/100 

YA=I Y-100*C 

J = ! 146097^=0 /4+! 1461*YA) /4+ ! 1 5 3=^ I M+2 ) /5 +D + 172 1 1 1 9 

RETURN 

END 



SUBROUTINE SHELL!N,X) 
DIMENSION X!l) 

M=N 

0101 M=M/2 

IF(M.EQ.O) RETURN 
K = N-M 
J = 1 
0111 I=J 
0121 IM=I+M 

IF!X! I ) .LE.Xn M) ) GO TO 0131 
S=X! I ) 

X! I ) =X! IM) 

X ! I M ) = S 
I =I-M 

IF! I .GE.l) GO TO 0121 
0131 J=J+1 

IF(J.LE.K) GO TO 0111 

GO TO 0101 

END 
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